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Abstract 


In this paper, we proposed a notion of belief universal gravitation (BUG) in the Dempster- 
Shafer (D-S) evidence theory, of which the notion of mass of a belief function is newly 
addressed using evidence quality coding (EQC) method. The proposed BUG aims to 
discuss the process of information fusion from the perspective of Newton’s mechanics, 
which may provide us a new insight to address the issues of D-S evidence theory. A key 
issue in D-S evidence theory, 1.e., conflict management, is solved better than previous 
methods using the proposed BUG. An application in fault diagnosis is used to illus- 
trate the effectiveness of the proposed BUG. Some further work is also summarized to 
present the potentials of the proposed BUG. 


Keywords: Information fusion, Dempster-Shafer evidence theory, Belief universal 
gravitation, Evidence quality coding algorithm, Conflict management, Fault diagnosis. 


1. Introduction 


Dempster-Shafer (D-S) evidence theory was first proposed by Dempster and then 
further developed by his student Shafer [1, 2]. In the 1970s and 1980s, D-S evidence 
theory was introduced into the field of artificial intelligence [3]. Like fuzzy sets [4, 5], 
rough sets [6, 7], Z-number [8, 9, 10], belief structures [11], D numbers [12, 13, 14], 
soft likelihood functions (SLF) [15, 16, 17], belief entropy [18, 19] and belief function 
[20] as an uncertainty reasoning method, D-S evidence theory provides a powerful 
tool for the representation and fusion of decision-level uncertainty information. Up to 
now, D-S evidence theory has been widely used in information fusion [21, 22], fault 
diagnosis [23], decision analysis [24, 25], and so on, due to its advantages in dealing 
with uncertain information. 

In particular, Dempster’s combination rule is the core part in evidence theory, which 
is used to fuse evidence information from multiple independent sources. However, the 
founder of fuzzy mathematics Zadeh found through a counter-example: when using 
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it to combine conflicting evidence, the result is contrary to human intuition [26]. For 
a deeper understanding, we briefly review this example. From Definition 5, it can be 
clearly observed that there is a clear conflict between the evidence mı and m2. The 
reason is as follows: 


e Proposition A is strongly supported by the mı, but completely denied by the mp. 
e Proposition C is strongly supported by the m, but completely denied by the mı. 


Then the Dempster’s combination rule was used to fuse the two pieces of evidence mı 
and m, and the fusion result was completely positive for proposition B, which was low 
in support of original evidence mı and m2. Obviously, this was not intuitive for human 
beings, and the fusion result was not convincing. 

Since then, more and more scholars have questioned the validity of Dempster’s 
combination rules. The research on the conflict between evidences has always been a 
hot topic in D-S evidence theory. Generally speaking, there are two ways to deal with 
conflict information. The first way is to modify the combination rule of the classical 
D-S evidence theory to adapt to the environment of high conflict, that is, to modify 
the rule. The second idea is to keep the combination rule of the classical evidence 
theory and preprocess the conflict evidence before fusion. In other words, this type of 
method is to modify the data model. Where, the viewpoint of modifying rule can be 
summarized as follows. Conflict management is a key problem to improve and develop 
evidential reasoning. The use of Dempster’s combination rule under high conflict of 
evidence will produce unreasonable conclusions, which is generated by the normaliza- 
tion step of the rule, so it is necessary to modify Dempster’s combination rule. The 
new combination rule needs to focus on how to redistribute conflicts. The representa- 
tive of this school is the ” Unified reliability function combination method” proposed 
by Lefevre et al. [27]. Other methods under this thinking include: Yager [28, 29], In- 
agaki [30] et al. introduced the method of average support for propositions, Sun et al. 
[31], Zhang’s method [32], the conditions of Dempster’s combination rules [33], the 
combination rule of minC proposed by Dinael et al. [34] refined allocation space and 
local conflicts and proposed the concept of potential conflicts, Smardndache et al. [35] 
proposed the conflict proportional allocation rule PCR3 and Ma et al. [36] recently 
proposed the flexible combination rule based on the complete conflict set, and so on. 
Unlike the first class of methods for modifying combination rule, the second class of 
methods is based on modifying the original evidence source. The view of this class 
of methods is that Dempster’s combination rules are not inherently wrong. When the 
evidence is highly conflicted, the conflict evidence should be preprocessed first, and 
then the Dempster’s combination rule should be used. Such methods are represent- 
ed by Haenni et al. [37]. Other typical methods include Murphy’s simple arithmetic 
average method [38], Deng et al.’s weighted average method [39], et al [40, 41]. For 
now, researchers lean toward Haenni’s point of view, and papers that modify data mod- 
els dominate. In particular, Liu [42] comprehensively considered the applicability of 
Dempster’s combination rule in the case of conflict, proposed the method of using bina- 
ry groups of k and dif BetP to describe the conflict, and proposed the proposal of using 
Demspster’s combination rule in the authoritative journal Artificial Intelligence. So 
far, research on conflicts in evidence theory continues [43, 44, 45]. 


Above, we briefly reviewed the development history and current status of D-S evi- 
dence. Since the proposed method is based on the theory of gravitation, it is necessary 
to introduce some basic knowledge related to the theory in the following paragraphs. 

The law of gravitation was first proposed by Newton in the book ” Mathematical 
Principles of Natural Philosophy” published in 1687 [46]. The law of gravitation be- 
longs to the laws of natural science. It shows that any two objects in nature are attracted 
to each other. The magnitude of gravity is proportional to the product of the masses 
of two objects and inversely proportional to the square of the distance between them 
[46]. Subsequently, the law of gravitation has been widely developed and applied in 
the field of natural science. In recent years, some scholars have proposed a population 
optimization algorithm based on the idea of the law of gravitation [47], while others 
have developed a gravitational search algorithm [48], and so on. This novel algorithm 
combining physical meaning has attracted the attention of many scholars immediately. 
Since then, more and more scholars have tried to apply gravitation to relevant scien- 
tific research fields. Therefore, studies based on gravity have been further developed 
[49, 50, 51]. 

Based on the above discussion and knowledge background, inspired by Newton’s 
law of gravitation, in this paper, for the first time, we creatively proposed the belief 
universal gravitation (BUG) in D-S evidence theory. First, we believed that the essence 
of evidence information fusion is affected by some potential force. If certain conditions 
are met, the evidence will be fused. Secondly, in order to quantify the degree of force 
between the evidences, we proposed the theory of belief universal gravitation (BUG). 
Among them, the BUG formula is the core part of the theory. The details about it are 
expressed as follows. 


e The evidence obtained by the sensor is abstracted into the quality of logical evi- 
dence. The proposed evidence quality coding (EQC) algorithm is used to obtain 
the quality of evidence for each independent source. 


e The evidence distance [52] is used to indicate the spatial distance between two 
pieces of evidence. 


e The evidence gravity parameter Ggr is used to distinguish the different discern- 
ment frameworks in the system where the same evidence gravity is located. 


In addition, by proof, the BUG formula satisfies some basic properties. The relation- 
ship between it and other variables is further explained through simulation experiment. 
Thirdly, based on the BUG theory, we modeled the essence of information fusion. In 
addition, combined with the clear physical meaning of the BUG formula, we studied 
its effect in conflict management. Finally, application illustrates the superiority and 
prospects of the proposed method. 

The organizational structure of this paper is as follows. Section 2 introduces some 
background knowledge required for this paper. In Section 3, the BUG theory is pro- 
posed, and some basic properties of the BUG formula are proved. In Section 4, a nu- 
merical example is used to model the essence of information fusion in evidence theory. 
Section 5 illustrates the nature of the proposed method through simulation experiment. 
In Section 6, the proposed BUG formula is used for the measurement of conflict, and 


the potential of the proposed method is illustrated by comparison. Section 7 shows the 
application value of the proposed method through the application of fault diagnosis. 
Section 8 summarizes and discusses the work done in this paper. 


2. Preliminaries 


In this section, we briefly review some basic concepts, including D-S evidence 
theory, evidence distance, and law of gravitation. 


2.1. Dempster-Shafer evidence theory 


2.1.1. Frame of discernment 
DEFINITION 1. Let © be a set of mutually exclusive and collectively exhaustive events 
defined by [1, 2] 
© = {01,02,--- , On} (1) 
where the set © is called the frame of discernment. 
The 2° is the © power set, which is expressed as 


2° = {0, {01}, +++ {Ov}, {01,02}, -++{01,02,--* ,Oj},--» O} (2) 


and Ù is an empty set. 
If A € 2°, A is called a hypothesis or proposition. 


2.1.2. Mass function 
DEFINITION 2. For a frame of discernment ©, a mass function is expressed as a mapping, 
i.e., from 2° to [0, 1], formally defined by [1, 2] 


m: 2° — [0,1] (3) 
which satisfies the following two attributes 


m(0)=0 and > m(6)=1 (4) 


OCO 


In D-S evidence theory, m is also called a Basic Probability Assignment (BPA). 
For exmaple, m(A) is BPA of A, which accurately reflects the extent to which A is 
supported. If m(A) > 0, A is a focal element of the mass function. 


2.1.3. Belief and plausibility functions 
DEFINITION 3. From this BPA, a belief function Bel and a plausibility function PI are 
defined, respectively, as [1, 2] 


Bel(A) = ` m(B) (5) 
BCA 
and 
PI(A) = 1 — Bel(A) = ` m(B) (6) 
AN B=0 


where A = ® — A, Bel : 2° > [0,1] and PI : 2° — [0, 1]. 


The relationship between P/ function and Bel function is shown in Fig.1. 


Min belief level of A Max belief level of A 


Bel(A) he level of ignorance in A Bel(A) 





Figure 1: The relationship between PI and Bel. 


In Fig.l, the quantity Bel(A) can be interpreted as a measure of one’s belief that 
hypothesis A is true. The plausibility P/(A) can be viewed as the total amount of be- 
lief that could be potentially placed in A. The [Bel(A), PI(A)] indicates the uncertain 
interval for A. 


2.1.4. Dempster’s combination rule 
The Dempster’s combination rule has been widely used to combine multiple inde- 
pendent evidence, and its definition is as follows. 


DEFINITION 4. Suppose the two evidence functions m; and m are on the discernmen- 
t frame ©, and then the Dempster’s combination rule can be defined as follows (® 
represents the orthogonal summation operation.) [1, 2] 


0 d= 
[m; ® m] (A) = : > Oe) (7) 
where the conflict coefficient k is defined as follows 
k= Š, mj(Ai)ma(A2) (8) 


Ai NA2=0 


Notably, Dempster’s combination rule is useful only under the condition that k < 1. 


2.2. Zadeh’s counter-example 


DEFINITION 5. Suppose m, and m are two BPAs defined on a frame of discernment 
© = {A, B, C} with [26] 


mı : m (A) = 0.99, mı (B) = 0.01, mı (C) =0 


m : m (A) = 0, m(B) = 0.01, m2(C) = 0.99 


Then using Dempster’s combination rule, the fusion result of proposition B is 
2 m(X)m(Y) 
XNY=B 


l- } m(X)m(Y) 
XNY=0 


Meim m} ({B}) = 


mı ({B}) m ({B}) 


1 = (m ({A})) m ({B} + m ((A1} m (AC) + mı (B) m (CH) 
0.01 x 0.01 


1 — (0.99 x 0.99 + 0.99 x 0.01 + 0.01 x 0.99) 
1 


2.3. Evidence distance 


Jousselme et al. [52] proposed a distance measure for belief functions, the evidence 
distance is defined as follows. 


DEFINITION 6. Let m; and m be two BPAs on the same discernment frame ©, and the 
distance between m, and mz is defined as follows [52] 








ALQA 
d(m,,m2) = |> ` E — m(A1))(mı (A2) — m(A2)) (9) 


0#A,CO Aj U A2] 
0#A,CO 


DEFINITION 7. A metric distance defined on the set Ñ is a function. Evidence distance 
d satisfies the following conditions [52] 


x> ÑK (10) 


(A, B) > dA, B) 
which meets the following requirements for any of A and B in Ñ: 
(1) Nonnegativity: d(A,B)> 0 
(2) Nondegeneracy: d(A,B)= 0& A=B 
(3) Symmetry: d(A,B)=d(B,A) 
(4) Triangle inequality: d(A,B)< d(A,C)+d(C,B) and any C € 6 


2.4. Law of gravitation 


DEFINITION 8. Let F be the magnitude of the gravitational attraction on the object, G 
be the universal gravitational constant, M,and M, be the masses of any two objects 
attracted to each other, and R be the distance between the two objects. And then, the 
formula for gravitation is defined as [46] 


MM) 
R2 





F=G (11) 


From Eq.(11), we can know that the law of gravitation is expressed as follows. 
Any two objects in nature are attracted to each other, and the magnitude of gravity is 
proportional to the mass product of the two objects, and inversely proportional to the 
square of their distance. 


3. Proposed the theory of belief universal gravitation 


In this section, inspired by Newton’s theory of gravitation, based on evidence dis- 
tance, we creatively present a novel BUG within the framework of evidence theory. 


3.1. Proposed the evidence quality coding algorithm 


In order to characterize the quality of evidence obtained from evidence sources, 
in this subsection, we propose the EQC algorithm to generate evidence quality. The 
introduction to the method is shown below. 


ASSUMPTION 3.1. Assume that mı and mz are two pieces of evidence defined on the same 
discernment frame © = {A, A2, A3}, and their BPAs are shown below 


mı : mı (A1) = 71, m(A2) = Jo, mı (A3) = D3 
m : m(Aı) = 91, m(A2) = p2, M2(A3) = 3 
where 7; + 0, 2 # 0, 03 # 0, yı #0, Yo # 0 and g3 # Q. 


Step 1: Assign each BPA a binary code. Each BPA was assigned a binary code based 
on the order of propositions in the discernment frame. The coding principle is as 
follows: for each BPA, the corresponding bit of the proposition in the discern- 
ment framework is labeled 1 and the rest is labeled 0. The binary code of m is 
shown below 


m,(A;) — 100, mı (A2) — 010, mı (43) — 001 
Similarly, the binary code of m is as is represented as follows 


M7(A1) => 100, m(A2) => 010, M7(A3) — 001 


Step 2: Convert the binary encoding of each BPA to decimal. So, m, is converted as 
follows 
m(A;) > 4,mı (A2) > 2,mı (A3) > 1 


Similarly, m is converted as follows 


m (A1) > 4,m(A2) > 2, m(A3) > 1 


Step 3: Generate the quality of each piece of evidence. The quality of the evidence m; 
is as follows 
E 4 x P +2 XxX +1 xD _ 4 x P +2 xX +1 xi 
a n 7 3 
Similarly, the quality of the evidence my is as follows 


M —4xg+2xXgotlxg 4x91 t+2Xg+1X ys 
————————————— G I 
n 3 


where n is the number of focus elements in BPAs. 


3.2. Proposed the formula for belief universal gravitation 


In the previous subsection, we described the process of generating evidence quality 
using the proposed EQC algorithm through a model. In this subsection, based on the 
evidence quality, and the evidence distance, we propose a new BUG formula. The 
formula is defined as follows. 


DEFINITION 9. Let my and m be two separate and different evidences on the same 
discernment frame ©. The BUG formula is defined as 


Mm, Mm 
F Bpa = Clg P @ (12) 
where the Grr is defined as 
Ger = 10°?! (13) 
with 
O0<d<l (14) 


In Eq.(12), Mm, and Mm, represent the quality of the evidence mı and m generated 
using the proposed EQC algorithm. d represents the distance between evidence from 
two independent sources. Ger is the evidence gravitation parameter used to distinguish 
different discernment frames. Where 6 is an adjustable parameter that satisfies the 
constraint of Eq.(14). It is used to dynamically adjust the size of the BUG in a system, 
so that the BUG can be observed more clearly. Simultaneously, it is important to note 
that the value of 6 in different BUG formulas must be consistent in a definite system 
(the definite system represents the same environment in which the belief gravitation is 
used, including the time environment and the space environment.). 

The physical meaning of the BUG formula is as follows. F gpa represents the grav- 
itation between two different evidences mı and mz. As can be seen from the Eq.(12), 
in a definite system, for two different evidences on the same discernment frame, the 
BUG 1s proportional to the product of their quality and inversely proportional to the 
square of the distance between them. Taking this one step further, the smaller the dif- 
ference between the two pieces of evidence, the greater the gravitation between them. 
Conversely, the greater the difference between the two pieces of evidence, or even the 
conflict, the smaller the gravitation between them. 


3.3. Basic properties of the proposed belief universal gravitation formula 


In this subsection, we discuss and demonstrate some of the basic properties that are 
satisfied by the BUG formula. 


ASSUMPTION 3.2. Assume mı and m are two different pieces of evidence on the same 
discernment frame ©. In a definite system, F gpa is the gravitational pull between 
evidence mı and m, Mm, and Mm, are the quality of the evidence m, and m2 generated 
using the proposed EQC algorithm, and Ger is the evidence gravitation parameter. 


The BUG formula has the following properties. 


(1) Non-negative 


Proof. First, from Eq.(9), it can be known that the evidence distance d > O, 
clearly, d? > 0. Secondly, through Eq.(13), we can know that G, > 0. Finally, 
it can be seen from the subsection 3.1 that the quality of evidence generated by 
EQC algorithm is non-negative, 1.e., Mm, and Mm, > 0. Based on the above and 
the Eq.(12), we can observe that F'gp, is non-negative. L] 


(ii) Symmetry 


Proof. In Definition 9, the value of 6 is determined because of the BUG in a sys- 
tem defined. Based on this, since mı and m are on the same discernment frame, 
the value of Ggr is determined. According to the evidence distance symmetry 
property (see Definition 7), 1.e., d(A, B) = d(B,A). In Eq.(12), Mm, Mn, is the 
product of the masses of evidence mı and m, hence, which also has symmetry. 
In general , the BUG formula has a symmetry, 1.e., Fgpa (A, B) = Fgpa (B, A). O 


(ii) Unbounded 


Proof. Since the BUG is in a definite system, the value of 6 is defined. As can 
be seen from the Eq.(12), the Fgp, is proportional to the product of the mass 
of evidence mı and m, and inversely proportional to the number of discernment 
frame |©| and the evidence distance d. It can be seen from Definitions 6 and 7 
that the evidence distance d is the scale to measure the difference of evidence, and 
its value is between 0 and 1. When the discernment frame is infinite, the value 
of evidence quality will increase exponentially due to the binary encoding used 
in EQC algorithm. Even if the evidence gravitation parameter Ggr weakens the 
value of BUG to some extent, the value of BUG will approach infinity with the 
infinite expansion of discernment frame elements. = 


4. Model the essence of information fusion 


In this section, through an example, we use BUG to model the essence of multi- 


sensor data information fusion in evidence theory. 


EXAMPLE 4.1. Suppose the discernment frame is © = {A,, A2, A3}. At T = Os, sensor 1, 
sensor 2, and sensor 3 obtained three pieces of evidence, namely mı, m2, m3, respec- 
tively. After that, evidence m4, ms were obtained by sensor 4 and sensor 5 at T = 5s, 
T = 10s respectively. Their BPAs are shown in Table 1. 


Table 1: The initial evidences obtained by the sensors. 


Time 

T = 0s 
T = 0s 
T = 0s 
PSs 
T = 10s 


Sensor 


Si: 
So 
S3 
S4 


S5 


mı(-) 


: M(-) 
: 3(-) 
: m4(-) 


: ms(-) 


m({A }) 


0.99 
0 
0.97 
0.95 


0.94 


m({A>}) 


0.01 


0.01 


0.03 


0.05 


0.06 


m({A3}) 


0 


0.99 


Obviously, the evidence m is a piece of conflict evidence because other evidence 
supports the proposition A; very high, while evidence m has no probability support 
for it at all. Under the proposed BUG, the fusion results of these five pieces of different 


evidence information are shown in Fig.2. 
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Final state: T=10s 


— 


\ 
N 


— 





Notes: <——> indicates a large BUG, which will cause the evidence to be fused. 


X indicates a conflict. 


Iya ct fam a sn meget ae ee neem’ le Vm cmt Tem stm A hw msi in > Rt sm Dame i Pam ey oar Nc at 2 ee 


: <-—-—> indicates a small BUG, which will not cause the evidence to be fused. 7 


Figure 2: Information fusion modeling of five different evidences based on BUG. 


In Fig.2, Fig.(a) shows five pieces of evidence information in the free state. It can be 
observed that when T = 10s, that is, all evidences are obtained, they are also scattered 
and have no relevance. In this case, information fusion is a simple combination of the 
five pieces of evidence information. However, the fusion results are not conducive to 
the final decision, because there is interference information m2. If the BUG system is 
considered, then at each phase of obtaining the evidence, the evidence will gradually 
merge the potentially fused evidence information under the force of the BUG. In this 
system, the quality of each piece of evidence information, the distance between them 
and the size of the BUG between them are different. 

In phase | (T = Os), for the three pieces of evidence obtained initially, it is clear that 
m2 is a piece of highly conflicting evidence information. Based on the proposed EQC 
algorithm, the quality of m2 is relatively small compared to the other two evidences; 
based on the evidence distance, the distance between m1 and m3 is also the closest, as 
shown in Fig.(b). After that, we assume that under the action of BUG, m1 and m3 are 
fused to form a new evidence information M’. Because of the increase of proposition 
reliability, its quality also becomes larger than the initial evidence m1 and m3. At this 
time, there is still a conflict between m2 and M’, and the BUG between them is not 
enough to make the two fuse with each other, as shown in Fig.(c). In phase 2 (T = 5s), 
when the evidence m4 is added to this system, the force relationship under the action 
of BUG is shown in Fig.(d). At this time, m2 is still a piece of conflicting evidence. 
The difference is that the distance between M’ and m4 increases due to the increased 
reliability of M’ for the proposition, but the BUG between them is still relatively large. 
After that, m4 and M’ are fused under the action of BUG to form new evidence M”, as 
shown in Fig.(e). Similarly, in phase 3 (T = 10s), Fig.(f) shows the force relationship 
diagram after the evidence m5 is added, and Fig.(g) shows the final state of the system 
under the action of BUG. On the whole, as the newly formed evidence M’, M” and 
M” support for propositions increase, their quality also increases. At the same time, 
the distances between | and m3, m4 and M’, and m5 and M” also gradually increase 
due to changes in reliability. 

Based on the concept of the proposed BUG, in this section we discuss the physical 
significance of evidence information fusion by multi-source sensors in evidence theory. 
This provides a new idea for information fusion. 


5. Simulation experiment 


In this section, the relationship between BUG and other parameter variables is fur- 
ther discussed through simulation experiments. 


EXAMPLE 5.1. Suppose there are 20 elements in the discernment frame ©, such as © = 
{1,2,3,--- ,20}. The BPAs of two different evidences mı and m are defined as follows: 


mı : mı(2, 3,4) = 0.05, m,(7) = 0.05, mı (©) = 0.1, mı (A) = 0.8 
My) : m (1,2, 3,4,5) = 1 


Here, assuming that A is not a set of constants, it can be changed by discernment 
frame. It starts at 1 and ends at 20, increasing by an order of magnitude each time. 
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Then, we record the data of the mass of evidence mı and m, the square of the distance 
of evidence d’, and the change of the BUG in Table 2. 


Table 2: The change data of the parameters in the BUG formula. 


A Ger Mn, Mm, d’ F BPA 
{1} 10~!° 136,908.775 1,015,808 0.6175 22.5226 
{1,2} 107° =189,337.575 1,015,808 0.4714 40.7982 
23] 107!° 215,551.975 1,015,808 0.3255 67.2748 
{1,2,3,4} 107° 228,659.175 1,015,808 0.1795 129.3848 
{1,2,3,4,5} 107!° 235,212.775 1,015,808 0.0175 1365.1000 
{1,2,--- , 6} 107° =. 238,489.575 1,015,808 0.1509 160.5915 
(ls 2.8457) 107° 240,127.975 1,015,808 0.2529 96.4475 
{1,2,--- ,8} 107° 240,947.175 1,015,808 0.3255 75.2007 
{1,2,--- ,9} 107° 241,356.775 1,015,808 0.3828 64.0488 
{1,2,---,10} 10°! 241,561.575 1,015,808 0.4295 57.1251 
{1,2,---,11} 10°! 241,663.975 1,015,808 0.4684 52.4087 
{1,2,---,12} 10°! 241,715.175 1,015,808 0.5015 48.9557 
{1,2,---,13} 10°! 241,740.775 1,015,808 0.5301 46.3212 
{1,2,---,14} 10° 241,753.575 1,015,808 0.5552 44.2339 
{1,2,---,15} 10°! 241,759.975 1,015,808 0.5774 42.5288 
{1,2,---,16} 107° 241,763.175 1,015,808 0.5975 41.1001 
{1,2,---,17} 10°'° 241,764.775 1,015,808 0.6156 39.8940 
{1,2,---,18} 10°! 241,765.575 1,015,808 0.6322 38.8475 
{1,2,---,19} 10° 241,765.975 1,015,808 0.6474 37.9356 
{1,2,---,20} 10° 241,766.175 1,015,808 0.6615 37.1283 


' In this system, set the value of the adjustable parameter 6 to 1/2. 


It can be seen from Table 2 that with the increase of the number of propositions 
in set A, the evidence quality of my), 1.e., the value of Mm, increases gradually. Since 
the number of propositions in m remains the same, hence, Mm, is a constant. Further 
more, the relationships between d and the size of set A, F gp, and the size of set A are 
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shown in Fig.3 and Fig.4, respectively. 
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Figure 3: The trend of d? with respect to |A|. 
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Figure 4: The trend of Fgpa with respect to JA]. 


From Fig.3, we can see that the square of the evidence distance d? has the oppo- 
site trend of the Fgpa with the increase of the number of propositions in set A. More 
specifically, when set A approaches set {1,2, 3,4, 5}, the value of d tend to be the low- 
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est. Conversely, when the value of set A deviates from the set {1, 2,3, 4,5}, the value 
of F gp, increases. Especially, d? shows a decreasing trend until the size of set A is 5. 
When the size of set A is 5, that is, A = {1,2,3,4,5}, d* reaches the minimum value. 
And then, as the size of set A increases, d? increases, but it’s less than 1. It can be seen 
from Definition 6 that evidence distance d is used to express the similarity between 
evidences. When the size of set A is 5, that is, A = {1,2,3,4,5}, mı and m have the 
greatest similarity, that is, d is the smallest, so d? is the smallest. After that, the number 
of propositions in set A increases, so does the difference between mı and mz, so d? 
increases. In addition, from Fig.4, we can also observe that when the value of set A 
is 5, that is, A = {1,2,3,4,5}, Fgpa reaches the maximum value. In combination with 
the physical significance of the BUG, the difference between evidence mı and m is the 
smallest, that is, evidence distance d is the smallest. In other words, the similarity be- 
tween evidence mı and mz is the highest. Therefore, the universal gravitation generated 
by mı and m is maximized. Thereafter, as the size of set A increases, the similarity 
between mı and m gradually decreases, and the value of the BUG also decreases. 

To sum up, through the experimental simulation, we further clarify the relationship 
between the BUG and its parameters. We can conclude that, in a system, and given 
the discernment frame, the BUG is inversely proportional to the square of the evidence 
distance. In other words, when the distance between two pieces of evidence is smaller, 
at this point, the BUG is greater, the more attractive the two pieces of evidence are to 
each other. 


6. Conflict management based on the proposed belief universal gravitation 


How to manage conflicts in D-S evidence theory is a difficult and challenging prob- 
lem. More crucially, in the management of conflict evidence, how to describe the sim- 
ilarity between the evidence is a key and important problem. In this section, based on 
the proposed BUG formula and its physical significance, here we creatively use it to 
describe the similarity between evidence. 

Before comparing the proposed approach with the others, let’s briefly review some 
of the work of some scholars on conflict management in the DST. In [52], Jousselme 
et al. proposed a distance measure for belief function. In [1], the conflict coefficient k 
was first used to represent the degree of conflict between evidences. However, in [42], 
Liu pointed out that the k doesn’t effectively measure the conflict between two pieces 
of evidence. A two-dimensional conflict model is then proposed, where the pignistic 
probability distance [53] and the conflict coefficient k are united to represent the degree 
of conflict. After that, Daniel [54] defined the plausibility conflict between BPAs. 
Lefevre and Elouedi [55] put forward a novel called the Combination With Adapted 
Conflict (CWAC) rule. This rule provided an adaptive weighting between Dempster’s 
rule and conjunctive rule based on Jousselme et al.’s evidence distance. Ma and An [56] 
considered a method to combine conflict evidence with a probabilistic dissimilarity 
measure. Song et al. [57] studied the correlation coefficient for the relativity between 
two BPAs and used it to measure the conflict. Since then, to measure the correlation 
degree between the two pieces of evidence, Jiang [58] considered a new correlation 
coefficient that takes into account both the non-intersection and the difference between 
the focus elements. Recently, a new belief entropy is presented as Deng entropy [59, 
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Conflict 


60]. Inspired by Deng entropy, Pan et al. [61] proposed a novel association coefficient 
of belief functions to measure conflicts between evidences. 

Here we use the Example 5.1 in Section 5 to compare the degree of conflict between 
the two groups of BPA with the above-mentioned different methods. The results are 
shown in Table 3 and Fig.5. 





The size of A 


Figure 5: Comparison of conflicts between different parameter values 


As shown in Fig.5, more specifically, the d and k, show a consistent trend for con- 
flict measurement. When the size of set A approaches 5, the measurement value of 
conflict is the lowest. When the size of set A deviates from 5, the conflicting measure- 
ments increase. DisSim is not monotonic until the size of set A is 5, so it is not an 
effective way to measure conflicts. CWAC is always kept at a low level, and it is insen- 
sitive to the conflict change. The classical conflict coefficient k is always maintained at 
0.5, which cannot distinguish the variation of evidence mı. The method of P/_C takes 
A value of 0 when set A is 5, which is unacceptable. The k, method has between 0 
and 5 elements in set A, which is more monotone. When the size of set A is 5, the 
collision value is smaller than the above mentioned methods, which is reasonable. As 


16 


A 

{1} 

{1,2} 

{1, 2,3} 

{1, 2,3, 4} 
{1, 2,3, 4, 5} 
{1,2,--- ,6} 
{1,2,---,7} 
{1,2,--- ,8} 
{1,2,--- ,9} 
{1,2,--- , 10} 
{1,2,---,11} 
{1,2,--- , 12} 
{(1,2,---,13} 
{1,2,--- , 14} 
{1,2,--- , 15} 
{1,2,--- , 16} 
{1,2,--- ,17} 
{1,2,-+- , 18} 
{1,2,-+- , 19} 
2, -- B20} 


k[1] 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 


0.05 


Table 3: Comparisons of conflict degree. 


d[52] CWAC[55] DisSim[56] PI-C[54] k,[58] 


0.7858 
0.6866 
0.5705 
0.4237 
0.1323 
0.3884 
0.5029 
0.5705 
0.6187 
0.6554 
0.6844 
0.7082 
0.7281 
0.7451 
0.7599 
0.7730 
0.7846 
0.7951 
0.8046 


0.8133 


0.0393 
0.0343 
0.0285 
0.0212 
0.0066 
0.0195 
0.0251 
0.0285 
0.0309 
0.0328 
0.0342 
0.0354 
0.0364 
0.0372 
0.0380 
0.0386 
0.0392 
0.0397 
0.0402 


0.0407 


0.3710 
0.4855 
0.3974 
0.3644 
0.3375 
0.4188 
0.6000 
0.6497 
0.6884 
0.7194 
0.7448 
0.7660 
0.7839 
0.7992 
0.8126 
0.8242 
0.8345 
0.8438 
0.8519 


0.8389 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.05 


0.7348 
0.5483 
0.3690 
0.1964 
0.0094 
0.1639 
0.2808 
0.3637 
0.4288 
0.4770 
0.5202 
0.5565 
0.5872 
0.6137 
0.6367 
0.6569 
0.6748 
0.6907 
0.7050 


0.7178 


' In this system, set the value of the adjustable parameter 6 to 1/2. 
2 In this example, in order to facilitate comparison with other methods, we reduce every 
BUG in the same proportion (1.e., each BUG is divided by 1,000) without affecting the 
properties of the BUG. 
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CppalOl] 


0.9555 
0.8950 
0.7787 
0.5292 
0.0302 
0.5133 
0.7511 
0.8691 
0.9278 
0.9571 
0.9717 
0.9790 
0.9827 
0.9845 
0.9855 
0.9860 
0.9863 
0.9866 
0.9869 
0.9875 


F Bpa 
0.0225 
0.0408 
0.0673 
0.1294 
1.3651 
0.1606 
0.0965 
0.0752 
0.0640 
0.0571 
0.0524 
0.0490 
0.0463 
0.0442 
0.0425 
0.0411 
0.0399 
0.0388 
0.0379 


0.0371 


an improved method based on Deng entropy, the change of cgpa is more significant 
than that of k,. On the whole, with the increase of the number of elements in set A, the 
proposed BUG changes in an opposite trend to other methods. 

According to the physical meaning of the BUG discussed in subsection 3.2, it de- 
scribes the magnitude of gravity between the two pieces of evidence. That is to say, 
the more similar the two evidences support certain propositions, the larger the BUG. 
Conversely, the greater the difference in support for certain propositions, the smaller 
the BUG. In this example, although d, k, and cgpa can also represent the change of 
conflict, it should be emphasized that the proposed method shows from a mechanical 
point of view that the BUG can also well reflect the difference of conflict evidence. 
Fig.5 graphically reflects this change trend. This also shows the potential of BUG for 
conflict measurement. 


7. Application in fault diagnosis 


In industrial production, how to effectively integrate multi-source evidence to give 
decision makers reasonable and accurate machine fault diagnosis results is a challeng- 
ing problem. In this section, we demonstrate the superiority of the proposed method by 
applying the BUG to actual fault diagnosis in [62]. A machine fault diagnosis model is 
shown in Fig.6. 





Data fusion 


Data fusion result 
set 





set 


Data processing 
terminal set 


Final decision 
results 


Figure 6: Machine fault diagnosis diagram. 
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7.1. Application background 

In the case of motor rotor fault diagnosis, three kinds of sensors are used to collect 
the feature information of acceleration, velocity and displacement. The data from these 
sensors are turned into BPAs, which are listed as follow Table 4. 


Table 4: The output of the multi-sensors. 


Mi A A2 A3 Aa 
my 0.06 0.68 0.02 0.04 
m 0.02 0 0.79 0.05 


m3 0.02 0.58 0.16 0.04 


As shown in Table 4, mı, m2, and m3 represent evidence from acceleration sensors, 
speed sensors, and displacement sensors, respectively. Motor rotors typically have four 
states, A; is normal operation, A; is unbalance, Az is misalignment, and Ay is pedestal 
looseness. Then a discernment frame, 1.e., © = {A;, A2, A3, A4} is established. From the 
information collected by the three acceleration sensors, it is clear that the information 
obtained by the sensor my is conflicting because it does not support Az misalignment, 
which is highly supported by the sensor mı and m3. 


7.2. Fault diagnosis based on the proposed method 

In order to fuse the information obtained by the three sensors, we adopt the method 
proposed in [39]. Note that instead of using evidence distance to represent the sim- 
ilarity between the evidence information obtained by the sensors, we use BUG. The 
method is briefly expressed as follows. 


Step 1: Construct a BUG matrix (BUGM), which reflects the degree of similarity be- 
tween the two evidence bodies. Suppose there are n pieces of evidence, the 


BUGM as follows: 
3 Fppa(m,m2) +++ Fppa(m,mj) +++  Fgpa(mi, mn) 
BUGM =| Fgpa(mi,m;) Fpgpa(mim) +- Fepa(m,m;) + Fepa(m,m,)| (15) 
Fppa(m,,) Fgpa(mn,m2) +++ Fgpa(mn,mj) >> y 


where Fgpa(mi, m;) represents the gravitation between the evidence m; and mj, 
and J means not defined. 


Step 2: Calculate the credibility of all evidence. The degree of support of the body of 
evidence is expressed as follows: 


n 


Sup(m)= J, Fapa(mi,my) (16) 


jiji 
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and the credibility degree Crd; of the body of evidence m; is defined as: 


S i 
Crd; = a (17) 
2 Sup(m) 
Step 3: Get the average weight value of all evidence. 
MAE(m) = X (Cra; x mi) (18) 


i=1 


Step 4: Use the Dempster’s combination rule Eq.(7) to fuse the average weight of the 
evidence n — | times. 


More specifically, in the concrete calculation process, using Eq.(17), the credibility 
of the three pieces of evidence obtained is expressed as Crd, = 0.4872, Crdz = 0.1088, 
Crd3 = 0.4041 respectively. Through Eq.(18), the average value of the three obtained 
evidences 1s shown in Table 5. 


Table 5: The averaged evidence. 


m(A)) m(A) m(A3) m(A4) m(©) 
0.0395 0.5657 0.1604 0.0410 0.1934 


Finally, the evidence was fused twice using the Dempster’s combination rule, and 
then the final results are shown in Table 6. 


Table 6: The final results after using Dempster’s combination rule. 


m(A) m(A2) m(A3) m(A4) m(@) 
0.0111 0.8860 0.0763 0.0116 0.0149 


7.3. Analysis and discussion 


The outcome of a fault diagnosis depends on the probability of supporting a fault 
event. In this paper, we set the threshold to 0.8 to make a decision. Table 7 shows the 
initial diagnostic results of three different types of sensors. Table 8 and Fig.7 show the 
fault diagnosis results after the fusion of fault evidences obtained from three sensors 
by different methods. 
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Table 7: The initial output of sensors and diagnosis results. 


Sensor type 
Acceleration 
Velocity 


Displacement 


A 
0.06 
0.02 


0.02 


Ay 


0.68 


0.58 


A3 


0.02 


0.79 


0.16 


A4 


0.04 


0.05 


0.04 


Diagnosis result 


uncertainty 


uncertainty 


uncertainty 


By comparing the fault diagnosis results in Table 7 and Table 8, we can clearly see 
the advantages of multi-source sensor data fusion. More specifically, in Table 7, if we 
consider only one feature of the machine’s working state, it is difficult to make a correct 
judgment decision. With a threshold of 0.8, these three single pieces of evidence show 
that they present uncertain answers before the evidence is combined. In addition, how 
to deal with the conflicts in the evidence combination is also very important. As shown 
in Fig.7, if we use Dempster, Murphy and Deng et al.’s method to combine the evidence 
obtained by the sensors, it is difficult for us to make a decision. Jiang’s method and 
X1iao’s method can determine that the failut of the equipment is unbalance, while the 


proposed method seems to make more accurate judgments. 


Method 
Dempster [1] 
Murphy [38] 

Deng et al. [39] 
Jiang [58] 
Xiao [63] 
Proposed 


' In this system, set the value of the adjustable parameter 6 to 1/4. 


Table 8: Comparisons of some existing methods. 


Ay 
0.0205 
0.0112 
0.0110 
0.0108 
0.0146 
0.0111 


A2 
0.5230 
0.6059 
0.7730 
0.8063 
0.8184 
0.8860 


A3 
0.3933 
0.3508 
0.1856 
0.1534 
0.1021 
0.0763 
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A4 
0.0309 
0.0153 
0.0139 
0.0134 
0.0187 
0.0116 


© 


0.0324 
0.0168 
0.0165 
0.0162 
0.0462 
0.0149 


Diagnosis result 


uncertainty 
uncertainty 
uncertainty 
unbalance 
unbalance 


unbalance 


0.9 


0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.1 i 
Qe e n a ae i EE ee es es A aa 
Al A2 A3 A4 © 


m Dempster m Murphy Dengetal. mjiang mXao m Proposed method 


Figure 7: The results of comparison with some existing methods. 


To sum up, the example verifies the validity and superiority of BUG in the case of 
conflict and shows the potential application prospect of proposed method. 


8. Conclusion and discussion 


In this paper, a new concept of BUG is proposed in the evidence theory system. 
Overall, the contributions of this paper are as follows. First of all, to the best our 
knowledge, this is the first time to introduce the law of gravity from natural science 
into evidence theory. Secondly, according to the knowledge framework of D-S evi- 
dence theory, the BUG formula is used to characterize the degree of gravity between 
evidences. Through the modeling of evidence information fusion, the rationality of the 
BUG in the fusion process is explained. Finally, by solving the hot issue of conflict 
measurement in D-S evidence theory, the potential application value of the proposed 
method is demonstrated. Furthermore, the application also illustrates the effectiveness 
of the proposed method. 

In addition, although we proposed the theory of BUG, there are still some short- 
comings that need to be further overcome. For instance, under the BUG, what is the 
critical point of evidence information fusion? In other words, what conditions are met 
before the evidence can be fused. As a method of measuring the amount of evidence 
information, is there a more perfect alternative method than EQC algorithm? how to 
consider more applications to increase the reliability of BUG, and consider its appli- 
cation to more practical problems, etc. In future work, we will consider introducing 
speed and acceleration to establish the BUG theory more completely. We also intend 
to apply the theory of BUG to evidence anti-monitoring, interference interception and 
measurement of transmission media. 
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Highlights 


1. Anew concept BUG is proposed in the D-S evidence theory. 
2. BUG provides a new perspective for information fusion in evidence theory. 


3. BUG formula is based on evidence distance, EQC algorithm, etc. 
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